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Inducing topological transitions by a time-periodic perturbation offers a route to controlling the 
properties of materials. Here we show that the adiabatic preparation of a non-trivial state involves 
a selective population of edge-states, due to exponentially-small gaps preventing adiabaticity. We il¬ 
lustrate this by studying graphene-like ribbons with hopping’s phases of slowly increasing amplitude, 
as, e.g., for a circularly polarized laser slowly turned-on. The induced currents have large periodic 
oscillations, but flow solely at the edges upon time-averaging, and can be controlled by focusing 
the laser on either edge. The bulk undergoes a non-equilibrium topological transition, as signaled 
by a local Hall conductivity, the Chern marker introduced by Bianco & Resta in equilibrium. The 
breakdown of this adiabatic picture in presence of intra-band resonances is discussed. 


PACS numbers: 73.22.Pr,73.43.-f,03.65.Vf,05.30.Rt 

Introduction. Recent experiments [T] have mapped the 
phase diagram of the Haldane model [5], a prototypical 
Chern insulator, by driving ultra-cold fermionic atoms 
in an optical honeycomb lattice periodically modulated 
in time. Since the early suggestion by Oka & Aoki [3] 
for a photovoltaic Hall effect in graphene, theoretical re¬ 
search aimed at studying topological transitions induced 
by an external periodic perturbation — the so-called Flo¬ 
quet topological insulators [4] — has been intense li- 
m. Preparing a topologically non-trivial Floquet insu¬ 
lator, out of a standard band insulator, requires passing 
through a phase transition point, where the bulk energy 
gap momentarily closes and edge-states start crossing it. 
It is usually assumed that this can be done by keeping 
the system arbitrarily close to its Floquet “ground state” 
(GS), a concept which we will clarify later on, provided 
the strength of the periodic perturbation is ramped-up 
in an adiabatic way [3 uni HB — realizing a general¬ 
ized form of quantum annealing (QA) [221 - EH] . In this 
paper we study the QA dynamics of the Haldane model 
across its topological transition, and that of periodically 
driven graphene-like ribbons, e.g., irradiated by a circu¬ 
larly polarized laser. We find that the topological tran¬ 
sition comes with an ingredient that makes it different 
from the Kibble-Zurek (KZ) paradigm [26l [27] describ¬ 
ing the crossing of ordinary critical points |25l|28l|29]: an 
exponentially small Landau-Zener (LZ) [3011311 avoided- 
crossing gap between edge states, which forbids edge- 
state electrons from adiabatically following the GS, no 
matter how slowly the critical point is crossed. The pecu¬ 
liarity of this QA dynamics is reflected in non-equilibrium 
currents flowing at the edges, which could be controlled, 
e.g. by a laser focusing on the edges. We also show how 
the change in the topology of the non-equilibrium state is 
effectively signaled by the dynamical counterpart of a lo¬ 
cal Ghern marker, introduced by Bianco & Resta [32] as 
an indicator of an equilibrium non-trivial bulk topology. 

Model and idea. Graphene-like systems display re¬ 


markable properties associated to the pseudo-spin-1/2 
A — B sublattice degree of freedom of the honeycomb 
lattice, with relativistic Dirac cones sitting at the two 
corners K± = ±|^) — a being the lattice constant 

— of the hexagonal Brillouin Zone (BZ), when inversion 
symmetry (IS) and time-reversal (TRS) are unbroken. A 
minimal single-orbital tight-binding model, allowing for 
a time-dependence in the nearest-neighbor (nn) hopping 
phase and in the on-site energy, is given by the following 
Hamiltonian (omitting spin indices): 

H(t) = + AAB{t)Y^i-iyclci , ( 1 ) 

(ij) t 

where cl creates a particle at site i, {ij) denotes sums 
over nn and (—1)* = - 1 - 1 /— 1 on A/B. Aab controls IS, 
opens up a trivial equilibrium gap at the Dirac points 
and is in principle controllable versus time in optical 
lattice experiments [T]. The phases ^ij{t) — generally 
breaking TRS — may result from a time-periodic modu¬ 
lation of the optical lattice in the neutral cold atoms ex¬ 
periments [ 1 ], or from the Peierls’ substitution minimal 
coupling of the electrons with the (classical) electromag¬ 
netic field of a laser ^ij(t) = ^ • A(x, t). With the 

latter realization in mind, we take the field monochro¬ 
matic and described by a vector potential A(x, t) = 
Ao(x, t) [xsin(a;t) -I- y sin(a;t — 1 ^)], where ip describes a 
general elliptical polarization of the laser and Ao(x, t) is 
a smooth function of space and time. 

For a circularly polarized {(p = 4zTTj2) spatially uni¬ 
form laser [31 jT] |T31 [H], (t) = A(t) sin(a;t -I- cfij), with 

(fij = (±|,=Ff,7’‘) along the 3 nn directions (di,d 2 ,d 3 ) 
connecting an A-site to its nn H-sites, and X{t) = 
d being the nn distance. If the frequency w 
is larger than the unperturbed bandwidth W = 6|ti|, 
and A(t) and AAB{t) are nearly constant during a pe¬ 
riod T = 2ttIu!, the resulting Floquet evolution operator 
17 (t, 0) = gg effective Floquet Hamiltonian 


2 


TL^ approximately given by a Haldane model Hh with 
flux (j)H = ± ^, the same on-site difference A.ab , and hop¬ 
pings renormalized by Bessel functions: ti —)■ tiJo(A), 
^2 = —■\/3[fi>/i(A)]^/(fiw) [?]• As the amplitude \{t) is 
slowly turned-on — and/or AAB{t) is slowly decreased to 
0 — we effectively drive the Haldane model HH{t) across 
its equilibrium critical point {Aab/ t 2 )c = 3-\/3 [ 2 ]. In 
what follows, we study zig-zag strips with open boundary 
conditions (OBC) and -A/r sites in the cc-direction, and pe¬ 
riodic BC (PBC) along y (see inset in Fig. For each 

of the Ny y-momenta A:, the single-particle Hamiltonian 
is a Nx X Nx matrix t) whose Schrddinger dynamics 
is numerically integrated with a 4‘^-order Runge-Kutta 
method, the initial state |d'(0)) being the Slater determi¬ 
nant GS of H{Q) at half filling [ 55 ] , 

As a warm-up, let us first consider a QA of the Hal¬ 
dane model. Its phase diagram 0, AAB/t2 vs (///, 
is shown in Fig. a, the shaded regions denoting the 
topologically non-trivial phases with Chern number C = 
±1. In the insets, we show three zig-zag spectra at 
(/ij = f • a trivial insulator with AabI^^ = 4-\/3, the 
critical point {AABlt 2 )c = 3-\/3, and the IS-symmetric 
point with Aab = 0. Edge states cross the bulk gap 
in the non-trivial phase; the crossing fc-point between 
the two branches moves from the bulk-projected Dirac 
point iF+ = 27r/(3a) towards Kf = n/a as AAB/t 2 de¬ 
creases from {AAB/t 2 )c to Aab = 0. Actually, this 
is an avoided-crossing LZ point, with an exponentially 
small gap ~ where is the strip width and 

^ the localization length of the edge states, separating 
the two quasi-degenerate edge states. Consider the evo¬ 
lution denoted by the arrow in Fig. [^a: AABlt 2 starts 
from 4-\/3 and ends in 0 in a time Tq*. In the initial 
part of the evolution, the Dirac-point (bulk) gap Ak^ 
closes at criticality as Akj- ^ IjL, resulting in a stan¬ 
dard KZ [26l [27] non-adiabatic excitation of electrons 
into the conduction band [25]. In 2d, the critical ex¬ 
ponents 12 = 1, z = 1 should lead to a residual en¬ 
ergy Eres{t) = (4'(f)|i7(t)|4'(f)) - Ags(t), Egs(t) being 
the instantaneous GS energy, scaling as eres = Ej^e^it = 
Tqj^)/L'^ ^ Tqj/. The bare data for £^5 (stars in Fig. 
c) depart from this KZ scaling, due to a mechanism 
of selective edge-state excitation which we now discuss. 
Consider the right-edge electron sitting immediately to 
the right of the LZ gap in Fig. b (here we have only 
Ny = 72 A:-points for clarity of illustration): as the LZ 
gap sweeps towards larger k, it will be unable to follow 
the ground state due to the exponentially-small LZ gap, 
and will remain in the right-edge band, see Fig. Sb', 
but excited, since the equilibrium lowest-energy state sits 
in the left-edge band. In essence: there cannot be any 
LZ tunneling across the opposite edges of the sample. 
Hence, left-edge states remain, one after the other, se¬ 
lectively unoccupied. If we separate the contributions 
due to bulk and edges, Eres = CbuikL^ -I- SedgeL, we find 
that Ebuik ~ filled squares in Fig. [^c, while Eedge 



FIG. 1: (Color online) (a): Phase diagram of the Haldane 
model, with three representative zig-zag strip spectra, along 
the path of the QA evolution (arrow), (b, b'): The mech¬ 
anism by which right-edge states get selectively populated 
as the exponentially-small LZ gap sweeps to larger values of 
k during the QA evolution. The (equilibrium) edge states 
shown refer to Aab It 2 = 2.5 (for b) and AAB/t 2 = 2.4 (for 
b'). Filled/empty circles denote occupied/empty edge states. 
The blue vertical arrow in b points to an occupied right-edge 
electron that remains occupied (red vertical arrow in b') after 
the Landau-Zener event, (c): Residual energy (stars), sepa¬ 
rated into bulk (squares) and edge (circles) contributions, vs 
the annealing time tqa . Here data with L = N^. = Ny = 6n 
from 18 to 102 were used to get Ebuik/edge for each tqa. 

(solid circles) slowly increases, approaching the value 
^elge = IkI^ [Ek,+ -Ek,-], where Efc,+/_ are the 
right/left final edge bands. Starting the QA evolution 
from negative AabI 1 ^ 2 ^ or having (pn = ~f > swaps the 
role of right and left, and of the two Dirac points. 

Adiabatic Floquet Results. We now return to our Flo- 
quet QA, Eq. 0> with a circularly polarized driving. 
We performed an “adiabatic” linear turning-on of A(f) = 
(f/TQA)Af for a time Tqa = u-qaT, followed by an evolu¬ 
tion with constant Af for Tf = UfT (with t = 27 r/a;). We 
considered both a time-independent Aab , and a Aab (t) 
switched-off to 0 during the annealing of X{t) (physi¬ 
cally possible in the cold-atom realization). In princi¬ 
ple, Aab = 0 in graphene, but the initial zig-zag edge 
states would be pathological symmetric/antisymmetric 
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FIG. 2: (Color online) (a): Final Floquet quasi-energy bands 
for a uniform driving with (f) = —7r/2, hui = 7\ti\ > W, 
Aab = 10“^|ti| (effectively representing graphene), and \{t) 
linearly ramped up to Af = 1 in tqa = lOOr, with r = 27r/a;. 
Here the topological transition occurs at Acr ~ 0. Valence 
states (in red) are filled, conduction states (in black) empty. 
Filled circles denote occupied edge states. Inset: the initial 
graphene spectrum at A(0) = 0. (a'): Time-averaged bond 
currents calculated with a periodic evolution at constant Af 
for T{ = lOOr. The inset shows the large intra-period oscil¬ 
lations. (b,b'): Same as a, a' for an inhomogeneous driving 
focused on the right edge {xc = L^) of width a = QAL^. In 
the inset, a sketch of the zig-zag strip. 


combinations of left-right wavefunctions [51]. Hence, to 
describe graphene we include a very small Aab = 0^ 
whose value is not crucial (only the sign matters); the 
results turn out to be identical to evolutions in which 
AAB{t) is switched-off to 0. A generalized adiabatic the¬ 
orem holds for Floquet systems jT] HD] HU 155] |55] : a Flo¬ 
quet state |^/>q(A(0))) evolves remaining close to the in¬ 
stantaneous Floquet state \ilJa{X{t))) for sufficiently slow 
variations of the driving amplitude A(t), compared to the 
gaps from neighboring states rninmezdca — e /3 + rnHuj\). 
For hu) > FF, the standard Floquet BZ [—huj/2,hu!/2] 
is such that the initial Slater determinant |'F(0)) coin¬ 
cides with the Floquet ground state |'I'fgs(A(0) = 0)): 
all negative quasi-energies, in [—hui/2,0), are occupied, 
the positive ones are empty. The only relevant gap for 
the adiabatic Floquet dynamics [201 El] is that at the 
Dirac points. A slow increase of A will reproduce the 
QA of the Haldane case, as we verified by monitoring 
the occupations nk^a of the instantaneous single-particle 
Floquet modes \4>k,a) [3S]: the state |5'(rQA)) after the 
annealing is “close” to |'I'FGs(A(rQA))), apart from bulk 
KZ excitations near the Dirac points, and the previously 
discussed selective excitation of edge states, see Fig. [^a. 

The dynamics of the Hall current is interesting. As 
customary, in a Laughlin cylinder geometry the total 
current in the y-direction is given by Jy = ^ 1 , 


where Ky = is related to the flux <1)^, pierc¬ 

ing the PBC-cylinder along the az-axis, and t/fg is the 
flux quantum m- With PBC along y, we can write 
H{t) = c\ kC., f^, where H(fc,t) is the 

fc-resolved strip Hamiltonian and c-^ creates an elec¬ 
tron of momentum k at site i along the zig-zag line 
sketched in Fig. Hb'. Hence, bond-resolved y-currents 
are J^,i+lit) = {^{t)\YT'‘ ^i^r+i{k,t) 
with J = ^g. The circles in Fig. j^a' denote 

the time-average of Jiy+i(t) during the constant-Af evo¬ 
lution, [Jj,i+i]av = /t'^a Ji,*+i(t). Currents 

are concentrated at the edges, but with large periodic 
oscillations, shown in the inset: the stroboscopic aver¬ 
ages are not representative of the true time-averages [55] . 
Notice that the edge currents have left/right symmetry, 
while one would naively expect currents only on the edge 
which is selectively occupied by the out-of-equilibrium 
dynamics (the right edge, for Fig. [^. This behavior 
originates from specific symmetries — previously noted 
in equilibrium for the Haldane model at (j)H = ±7r/2 
[I1[IH] — whereby the GS value of Jy is zero everywhere 
due to an exact compensation between currents due to 
edge states and edge-current contributions due to bulk 
states. Since the out-of-equilibrium dynamics brings a 
lack of current-carrying edge states (at left, in Fig. 
the corresponding bulk contribution is uncompensated 
and gives rise to a left-flowing current of the same sign 
and amplitude as that at the right-edge. This left/right 
symmetry can be removed by a space inhomogeneity of 
the perturbation, e.g., a laser focused off-center. We find 
here expedient to retain translational invariance along y, 
assuming a y-independent Gaussian modulation ampli¬ 
tude Ag(x, t) = Ao{t) ^ where Xc is the focus 

center, and a the beam width. All the previous results 
remain valid for a central focusing, Xc = Ta;/2, provided 
a is not too small (cr > QAL^). But the interesting new 
feature is the ability to control the edge current to flow on 
either edge of the sample by moving the focus off-center. 
Fig.[2j-b illustrates the final Floquet quasi-energy bands 
when the laser is focused on the right edge {xc = L^), 
with a = OAL^. Notice that only the irradiated right- 
edge states show a A:-dispersion: unirradiated left-edge 
states stay flat and carry no current. Non-equilibrium 
currents flow only at the irradiated edge, see Fig. §b'. 

It is interesting to address the issue of an “indicator” of 
non-trivial topology in a non-equilibrium translationally- 
non-invariant setting. For translational invariant systems 
(with PBG), it was shown that the usual Ghern number C 
is conserved during a unitary evolution [11 [TH]: it does 
not “signal” the topological transition. However, since 
for fvjj > W the non-trivial final bulk states are adia- 
batically populated in a controlled way, we would expect 
to be able to “see” the topological transition by look¬ 
ing only at the “bulk” of the sample. We find that the 
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FIG. 3: (Color online) The bulk average C'buik(i) of the local 
Chern marker C(r,t) for a uniform driving with (p = —7r/2, 
hu) = 7|ti| > W, Aab = 0.1|ti|, and A(t) linearly ramped up 
to Af = 1 in tqa = 300t, followed by a constant-Af evolution 
for Tf = 220r, with r = 27r/a;. The topological transition 
occurs at Acr ~ 0.57. Here L = = Ny = 48, and we 

average on a central square of size 12 x 12. The horizontal line 
at « 0.96 is the time-average COuik, calculated from t = tqa 
to t = Tqa + Tf. The upper inset shows the saturation of 
Cbuik(-f '5 t"qa —>■ cjo) to a limiting value that, see lower inset 
(where we fit points with standard power-law corrections 1 -I- 
ci/L + goes to 1 for L —>• oo. 


local Chern marker C(r) introduced in Ref. |32|at equilib¬ 
rium, essentially a local measure of the Hall conductivity, 
works also in our non-equilibrium context: it signals if 
the sample bulk is locally a topologically non-trivial in¬ 
sulator, C(r) ^ ±1. C(r) can be expressed as a physically 
appealing commutator of position operators [32l |33] : 


C{r,t) = - 27 rf(r| [S:r{t),yvit)] k) ■ (2) 

Here x-p = VxV and yp = VyV are position op¬ 
erators projected on the occupied states, V{t) being 
the projector on the time-evolved Slater determinant 
Fig. I illustrates the dynamics of C(r, t), av¬ 
eraged over a central “bulk” portion of the sample, 
Cbuikk) = ^buikl^rebuik^k’*) a® the system evolves 
from a trivial insulator at A = 0, towards the non-trivial 
point with Af = 1. Upon time-averaging the oscilla¬ 
tions after Tqa [33], we obtain a quantity C§)([j.(L, Tqa) = 
C'buik(i) approaching the correct integer 
value 1 as Tqa —>■ oo and L ^ oo (see insets in Fig. |^. 

When huj is smaZ/er than the bandwidth W, intra-band 
resonances between valence states at €a and conduction 
ones at « Cq, + fiui change the physics completely, 
breaking the Floquet adiabatic picture. The whole crux 
is, in essence, that the starting point at A(0) =0 — 
the usual Slater determinant |'I>(0)) — does not coin¬ 
cide with the Floquet GS, |'I'(0)) ^ |'I>fgs(A(0) = 0)), 
for any choice of the Floquet BZ: upon folding the orig¬ 
inal [—lU/2, kF/2] spectrum in e.g., [—/ia;/2,/ia;/2], we 
see that the lower Floquet quasi-energy band is only 
partially filled, with empty conduction-band-originated 



FIG. 4: (Color online) Final Floquet quasi-energies and 
occupations nk,a (dot size proportional to rifc^^) for 4> = —7r/2, 
hijj — 4|fi| < W, Aab = 0.1|ti| (constant), and A(t) linearly 
ramped up to Af = 1 in tqa = 300r. 


states intermixed with filled valence-band ones. These 
filled-empty pairs with small Floquet gaps HOI m] 
rninmezdca — e /3 + mlkul) lead to a proliferation of bulk 
LZ events as A(t) is ramped up, with a complex redis¬ 
tribution of electronic occupations of the states: the fi¬ 
nal Floquet GS |4 'fgs(''^(''‘qa)) is nefer reached, even for 
Tqa —t OO. The analysis of these issues will be the sub¬ 
ject of a future work jSSj. Fig. shows the final Floquet 
quasi-energy bands, with the corresponding population 
indicated by a variable-size dot, for Iku = 4|ti| <W and 
X(t) linearly ramped up to Af = 1 in Tqa = 300r. The fi¬ 
nal state reached is a non-equilibrium metal, rather than 
an insulator. This aspect is important in devising pump- 
probe photoemission experiments in graphene |16j . 

Conclusions. We found a non-equilibrium mecha¬ 
nism which selectively populates edge states when per¬ 
forming an adiabatic switching-on of a periodic per¬ 
turbation towards a topologically non-trivial insulating 
phase. It is different from the “topological blocking” of 
Ref. |3H], which works with PBCs and when driving sys¬ 
tems with symmetry-protected subspaces from the topo¬ 
logically non-trivial phase to the trivial one. The mech¬ 
anism we illustrated requires edge states (hence OBG) 
whose electronic occupation is unable to follow instan¬ 
taneous equilibrium as they become topologically non¬ 
trivial and cross the bulk gap. It is general enough, 
and is at the root of the deviations from KZ scaling in 
Id topological transitions, as seen in |40l |4T]. In the 
present 2d context, it adds flexibility to the control of 
the edge currents flowing at the boundaries of the sam¬ 
ple, including the ability to have currents flowing only at 
one edge, by appropriate focusing of the ac field. Finally, 
we have shown that for huj < W intra-band resonances 
ruin the adiabatic picture and the resulting state is a non¬ 
equilibrium metal. Our findings should be amenable to 
experimental tests both with ultra-cold atoms in optical 
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lattices [T], and with laser irradiated electronic systems. 
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